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Abstract 

Harrell's modified perturbation theory [Ann. Phys. 105, 379-406 (1977)] is applied 
and extended to obtain non-power perturbation expansions for a class of singular 
Hamiltonians H = --£^ + x 2 + ^ + ^, (A > 0,a > 2), known as generalized spiked 
harmonic oscillators. The perturbation expansions developed here are valid for small 
values of the coupling A > 0, and they extend the results which Harrell obtained for 
the spiked harmonic oscillator A = . Formulas for the the excited-states are also 
developed. 
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1. Introduction 

This is a detailed extension of Harrell's modified perturbation theory 1 for the 
class of singular potentials 

H=-^+x* + ± ( A >0,a>2), (1.1) 

defined on suitable domains in the Hilbert space L 2 (0, oo) with solutions satisfying 
Dirichlet boundary conditions. By 'singular' we mean that the familiar Rayleigh- 
Schrodinger series either do not exist or do not converge. The present work, motivated 
by Harrell 1 and Greenlee 2 , studies a perturbative and variational analysis of the 
eigenvalues and eigenfunctions for the family of singular Hamiltonians 

H = H Q + \V = -^+x 2 + ^ + ± (A>0) (1.2) 
dx z x z x a 

known as generalized spiked harmonic oscillator Hamiltonian 3-10 . The extension lies 
in considering A to range over all non-negative real numbers instead of non-negative 
integers of the type 1(1 + 1). The main results are the extensions of Harrell's pertur- 
bative expansions 1 for the ground-state eigenvalues of the spiked harmonic oscillator 
Hamiltonian A = . In his elegant investigation, Harrell mentioned briefly the possibil- 
ity of extending his theory to the case of A = 1(1+1) , where I is the angular-momentum 
quantum number; however, his results mostly concern perturbation expansions for 
ground-state energies of the spiked harmonic oscillator Hamiltonian (1.1). There are 
two principal reasons for this choice: 1. The interesting Klauder phenomenon 11-13 
occurs only in the case A = , to the effect that, for sufficiently singular potentials, 
the perturbation term V cannot be smoothly turned off ( A — > ) in the Hamiltonian 
H = H +\V to restore the free Hamiltonian H ; 2. Rayleigh-Schrodinger perturbation 
series diverge at some finite order whenever a > 2 . 

Klauder's phenomenon doesn't occur 14-16 if A > 0. This is the case, for example, 
in N dimensions with A = (I + \(N - l))(l + \(N - 3)) and I > 0, or with I = and 
N ^ 1 or 3. In such cases the domain of the Hamiltonian H is stable under the limit 
A — > 0. However, a perturbative analysis for solutions that vanish at the origin is still 
interesting because of the divergence of the Rayleigh-Schrodinger series at some finite 
order for any a > 2 . We are able to conclude in the present article that the Rayleigh- 
Schrodinger series will breakdown at the order n > 2^(7 - 1) for a > 2 where v = 
and 7 = 1 + + AA . For example, a > 27 causes the perturbation series to diverge 
at the first order; for a > 7 + 1 the second-order perturbation will diverge, etc. These 
results and some others concerning the convergence of Rayleigh-Schrodinger series 
which rely heavily of the application of Kato's criterion 17-19 will be the subject of an 
independent investigation. In the present paper, we concentrate on the development 
of non-power perturbation expansions for the Hamiltonian (1.2). 

Detwiler and Klauder 20 , in their variational study of the spiked harmonic oscil- 
lator Hamiltonian (1.1), have shown that for 2 < a < 3 the eigenvalues are given by 
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asymptotic series to first-order so long as A > . But for a > 3 , the ground-state 
eigenvalues are given by 

E (X) = 3 + k\» + o{\ v ), 

and, for a = 3, by 

So (A) = 3 + fc'Alog(A) + 0(A), 

where k and k' are to be determined by variational means 20 . Harrell, soon af- 
terwards, modified the Rayleigh-Schrodinger series by utilizing the standard WKB- 
approximation technique for the lowest few orders. This proved to be quite successful, 
and he continued to developed a special perturbation theory, now known as 'singu- 
lar perturbation theory', and obtained thereby the first few terms of the perturbed 
A -expansion for different values of a . This turned out to be a non-power series ex- 
pansion and in fact was of exactly the same order as that of Detwiler and Klauder 20 . 
More specifically, Harrell 1 showed that the asymptotic series for the ground-state 
eigenvalues of the Hamiltonian (1.1) are explicitly given, for v = , by 

for a > 4 : 

W) - 3+ v^VTV) ( } ' 

for 3 < a < 4 : 

EoW = 3 + ^-"h" - jgSt> A + O(A-). 

for a = 3 : 

E (\) = 3 - -^=Alog(A) - ^A + 0(A 2 log 2 (A)) (c = .57721 56649 . . . Euler's constant). 

V7T y/TT 

For | < a < 3 : 

F m - 3 4- l^£(i^) r , g£(gg) 

The paper is organized as follows. In Section 2, we briefly review the regular per- 
turbation expansions for the Hamiltonian (1.2) and we identify the conditions un- 
der which the first- and the second-order corrections of Rayleigh-Schrodinger series 
exists. In Section 3, the main theorem used for the development of non-power per- 
turbation expansions will be introduced and proved. In Section 4, we introduced a 
suitable trial wave function. In Sections 5 and 6, we extend Harrell's theory to treat 
the generalized spiked harmonic oscillator Hamiltonians (1.2) for the cases where the 
Rayleigh-Schrodinger series fails and thereby we show that asymptotic series for the 
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eigenvalues of the Hamiltonian (1.2) are explicitly given by 
For a > 2(2 7 -l): 

^)-^ ^ T "rg^- a o)' 1) ' ^ , ' +0( ^ 1 '> 

For 2 7 < a < 2(27-1): 

£ (A) = 2 
For a = 27: 



p m - o,- 1 4(7 ~ 1 )^ (7 " 1)r ( 1 - Mr - 1)) aM7 -d 2i/r( 7 -£) n r x4^-i)x 
W) - 27 + r( 7 )r(i + 2K7-D) A "(1-2.(7-1)) r(7) A + ° (A } 



-c(l+7) + 21og(2(7-l)) 



(7 - i)r( 7 ) 



A + 0(A 2 log 2 (A)) 



^o(A) = 2 7 - 7 ---^ yf7 - y Alog ( A)- 
For 7 + 1 < a < 27: 

F (X) -2,+ 2^(7-Dr(l- 2^(7-1)) ^-!) , 2^(7-^) 2 

So(A) - 27+ ,r(7)r(2,(7-D) + (2K7-D-D r( 7 ) A + ° (A } ' 

where z/ = and 7 = \ + \\f\ + AA . The asymptotic expansions for the case a < 7 + 1 
are discussed in Section 7, along with some other cases. The connection with the region 
< a < 5/2 , overlooked by Harrell, is also investigated. In Section 8, the extension of 
the perturbation expansions developed in Sections 5, 6, and 7 to the excited states is 
discussed and some explicit formulas are derived. 

2. Asymptotic perturbation expansions 

It is known that although many perturbation expansions diverge, they may ac- 
tually be asymptotic expansions whose first few terms can yield good approxima- 
tions. The class of singular Hamiltonian (1.2) affords interesting examples of this 
phenomenon. Indeed, by regarding the Gol'dman and Krivchenkov Hamiltonian 7 
Ha = + x 2 + £ , which admits the exact solutions 



M*) = (-1)" J^-^-h-^^i-n^x 2 ) (2.1) 



with exact eigenenergies 



E n = An + 2 1 , n = 0,l,2,..., 7 = 1 + -V1+4A, (2.2) 

as the unperturbed part, and V(x) = x~ a as the perturbation potential, the first-order 
correction of the Rayleigh-Schrodinger series for the Hamiltonian (1.2) exists only for 
a < 27 , while the second-order correction will required a < 7 + 1 . The first condition 
a < 27 follows from E x = (ip ,x~ a ip ) , while the second condition a < 7 + 1 follows 5 
from 

" \fa,x- a 1>i)\ 2 



E 2 =Y 

^ Ei-E 

Under these conditions the perturbation expansions for the ground-state eigenvalues 
up to second-order 5 reads, for small values of A , 
E(X, a)= E + Si A + E 2 X 2 + ... 

- * + - ^(1, 1, 1 + f , 1 + \ ; 2,2,, + 1, 1) + . . . . ,2 ' 3> 
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In the next sections we develop non-power perturbation expansions for the cases where 
the regular Rayleigh-Schrodinger series fails to exist; namely a > 2-f and a > 7 + 1 . 

Before we proceed we should note that the functions 1F1 and 4 F 3 , mentioned 
above, are special cases of the generalized hypergeometric function 21 

p 

00 n (oti)k zk 

p F q (ai,a 2 , . . .,a p ;P 1 ,p 2 , . . . , q ; z) = ^ % ~T 77' ^ 2A " ) 

*=° n(&)fe 

where p and q are non-negative integers, and none of the fy, (j = l,2,...,q) is 
equal to zero or to a negative integer. If the series does not terminate (that is to 
say, none of the a; , i = 1,2, . . . ,p , is a negative integer), then the series, in the case 
p = q + 1 , converges or diverges accordingly as \z\ < 1 or \z\ > 1 . For z = 1 , the series 

is convergent provided J2 Pj ~ J2 a i > °- 

j=i »=i 

Here (a)„ , the shifted factorial (or Pochhammer symbol), is defined by 

(a) = l, (a)„ = a(a + l)(a + 2)...(a + n-l), for n = 1, 2, 3, . . . , (2.5) 

and may be expressed in terms of the Gamma function by (a) k = T(a + k)/T(a), when 
a is not a negative integer -m, and, in these exceptional cases, (-m)k = if k > m 
and otherwise (-m) fe = (-l) fe m!/(m - fc)!. 

3. The Main theorem 

It is clear that the perturbation approach mentioned in Sec. 2 cannot apply if 
a > 27 , since it is clear in this case that the first-order perturbation correction di- 
verges. We construct a modified perturbation series for the operators in this region 
by considering the perturbation theory of families of self-adjoint operators by an ap- 
plication of the variational method. This is done via Kato's generalization 22-23 of 
Temple's inequality 24-26 , which can understood from the following discussion. The 
derivation of bounds on the eigenvalues for self-adjoint operators usually starts from 
a consideration of the positive definite function given by 

(ji, M ) - {[H - £]</>, [H - e]4>) - (H4>, H<f>) - (0, H<t>) 2 + (e - (<f>, H<f>)) 2 > 0, (3.1) 

where fi is a function of 4> and e , i.e. \i = (i,(<f>, e) , H is the operator in question, e 
is a positive parameter, and 4> is a suitably chosen normalized trial function. If we 
expand the normalized function 4> in terms of the complete set of eigenfunctions {<?!>„} 
of H with eigenvalues E n (X) , <j> = J2 n a n<t>n , a n = {<j>, <j> n ) , (<j>, <f>) = 1 = |a„| 2 , we can 
express the positive definite function in (3.1) as 



( M , M )=]TK| 2 (£„(A)-e) 2 >0. 

n 
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Let us assume that we have picked the value of e for the closest approach to the ith 
eigenvalue E t , i.e. 

(p, M ) = £ K\ 2 (E n (X) - e) 2 > (Ei(X) - ef > 0. (3.2) 

n 

By combining (3.1) and (3.2), it can be easily seen that 

e - VII W1P#+(HP# < Ei(X) < e + y/\\ H<l>\\* -(<!>, H<l>)* + (e- & Ht))* 

(3.3) 

Now, by setting 26 

6 + [(H<t>, H4>) - (0, Hctif + (e - (<f>, H^)) 2 ] 1 ' 2 = Ef +1 (A), (3.4) 

where E^ +1 {\) is a lower bound estimate of E i+1 (X) , we can show that equation (3.4) 
possesses the solution 



1 

e= 2 



Ef +1 {\) -{<!>, H<f>) 



(3.5) 



provided (<p,H<p) < E[ +1 (X) . Substituting (3.5) into the lower bound expression in (3.3) 
yields the Kato- Temple expression for the lower bound: 

Ei{\) > (<p,H<p) - EUx) _^ m • (3-6) 

Similarly, setting 

e - [{H<t>, H4>) - (0, Hcji) 2 + (e - (0, tf0)) 2 ] 1/2 = iCi(A), (3.7) 

where £'^ 1 (A) is an upper bound estimate to the next lowest eigenvalue to Ei(X), 
yields 

for EjLiiX) < {<j>,H4>) . We let r\ = (<j>,H(j>), and the residual norm e =|| (H — ri)<j> | (hence 
e 2 =|| H<t> || 2 -r/ 2 ), and e 2 < (Ef- +1 (\) - r?)(r? - £^(A)), which follows by means of the 
inequalities r? - e 2 /(E l L +1 (X) - rj) > EV_ X or r) + e 2 /(ry - EV_ X {X)) < Ef +1 {X). This indeed 
ensures that the open interval (^ 1 (A),£'^ 1 (A)) contains a single isolated eigenvalue 
and no other piece of the spectrum. Then it follows from (3.6) and (3.8) that 



e 2 , , e 2 



This formula is symmetric with respect to upper and lower bound, as we might expect. 
It should be noted that (3.9) gives E^X) within the error bound of the order e 2 : this 
is very small if e is small, i.e. if is a good approximate eigenfunction. Indeed, (3.9) 
implies 

\Ei(X)-r]\<j, (3.10) 
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where g = min{r/-£'^ 1 (A), E^ +1 (X)-rj} . Therefore, the error in r\ depends on the residual 
norm squared, i.e. on e 2 , and on the gap g > for the eigenvalue is isolated 1 . If (3.9) 
is applied to the operator H + XV - , where Ef is a variational estimate for the 
i th eigenvalue of H Q + XV , there results 1 

Theorem 1: If $ is normalized trial function for the self-adjoint operator H = H Q + 
XV , where H and V are self-adjoint and Ef is an isolated, nondegenerate stable 
eigenvalue of H , and Ef is a continuous function such that (<j>, [H + XV - E*]<j>) — > 
as A — > , and 

|| [H + XV- Efo ||- o((<f>, [H + XV- E^) 1 ' 2 ), (3.11) 
then the eigenvalue of H + XV which converges to Ef satisfies 

Ei(X) = (<f>, [H + XV]4>) + 0(|| [H + XV- E*]<f> \\ 2 ). (3.12) 



Proof: To keep the notation simple, let us refer to EY_ 1 (X) and E^ +1 (X) in (3.9) 
by a and , respectively. Then, from the previous discussion, we have, for r\ = 
((f), [H + XV - Ef]4>), that a < r\ < (3, and e 2 < ((3 - r))(r) - a) . Further, by applying 
the Kato- Temple inequality (3.9) to the Hamiltonian H + XV - E* , we obtain, for 
normalized 4> , 

2 2 

V--^<E t (X)-E^<r,+ ^ ' 



a — n v i — / ' 

p — 7] T] — a 

where e 2 =\\ [H + XV - E*}4> || 2 -(</>, [H + XV '- E^]4>) 2 . If we divide by rj , we obtain after 
some simplifications 

\\[H +\V-E*]<j>\\ Y ( \\[H + XV-E}]4,\\ X 



' ) 71 Ej (A) — Ef 1 



< — ^ !__1<_^ L . (3.13) 



However, since is assumed to be normalized, 

E l (X)-El _ = E^-^Hq + XV]^^) = g i (A)-([ff + Ay]<M) / || [Hp + XV - E$}4> 
V V II [H + XV - E^}cj> || 2 ^ r?V2 

From (3.11), we have for A sufficiently small, 

|| [H + XV- E*]<f> 



v l/2 



< 1. 



Thus after dividing (3.13) by ^&«^Mj , we have 
Ei(X)-{[H + XV]4,,4>) 



[Ho + XV-E^\\ 2 



V I r 1 1 

- 1 1 ; — ( x max{- , } 



\\[H +\V-E^]4,\\ 



^1/2 



' (3 — rj 1 rj — a J 



< max{ — — — , — - — } < C (constant), 
p — 77 rj - a 
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which leads to 

Ei{\) = ([H + W]4>,4>) + 0(|| [H + XV- E*]4> f), 

as required. □ 

4. Trial wave function and solution to a differential equation 

In this section we shall introduce a suitable trial function in order to obtain eigenvalue 
perturbation corrections by means of Theorem I. For singular Hamiltonians of type 
(1.2), the trial functions are characterized by wave functions with non-integer expo- 
nent. This indeed characterizes 27-34 almost all trial functions which have been used 
previously to study this type of singular Hamiltonian (1.1) and (1.2). Furthermore, 
the trial functions have to satisfy the physical initial conditions of the problem. In 
the classical Rayleigh-Schrodinger perturbation theory, the lowest-order trial function 
for a given eigenvalue is chosen to be the unperturbed eigenfunction, i.e. the exact 
solutions of the unperturbed Hamiltonian. This is no longer a good choice for the 
perturbation XV in (1.2) with a> 2j , for 



oc 

/ 



x 2 '<- a - l e- x2 dx w t - a+2 '< 



which approaches oo as e goes to zero. Intuitively, it seems that if the unperturbed 
eigenfunction was modified slightly near the singular point, so that the expectation 
value of singular term V was no longer infinite, it would become a reasonable trial 
function to use to estimate the perturbed eigenvalue. This was the basic idea of the 
trial wave function used by Detwiler et al to study the Hamiltonian (1.1) and it 
was employed later by Harrell 1 . Using the notation of Harrell, we start with the 
(un-normalized) trial wavefunction 

i>(x;\) = W a (x;\)i>i(x), (4-1) 

where in(x) is given by (2.1) and W a (x;X) is to be determined. It should be noted 
that far away from the singularity, we expect ip(x;X) ~ ipi(x) for large x, since (1.2) 
behaves as radial harmonic oscillator Hamiltonian for large x , which, in turn, implies 
lim W a (x;X) = 1. Further, for an arbitrary singular point x , not necessarily at the 
origin, tp(x ;X) = , an idea that was borrowed from hard-core problems in quantum 
mechanics 20 : this forces W a (x ; X) = 0; therefore, we must also have lim W a (x; A) = . 
Using the trial function (4.1), the differential operator (1.2) leads to 



[H + XV- Ei]ip{x; X) = 



^w { x-,x)_MX)d +xvw • 

dx z dx dx 



^(x), (4.2) 



where Ei = Ef is the variational estimate of H . It is clear from (2.1) that « 
l^-ipiix) near the origin. Therefore, we may choose W a {x;X) in (4.2) such that 

* W f?X> + 2 fr-i)^»(*;A) _ WWa{x , A) = o (4 .3) 
dx z x dx 
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and must satisfies the initial conditions 

lim W a (x; A) = 0, and lim W a (x; A) = 1. 

X^O X^QO 

Eq.(4.3) allows us to write Eq.(4.2) as 

t dW a (x;X) 



[H + XV - E^j(x; X) = 2- 



dx 



x dx 



ipi{x) 



(4.4) 



(4.5) 



To solve (4.3) explicitly, we notice first that the parameter A can be removed from 
the equation by a change of variable z = X~ v x where v is to be determine shortly. A 
straightforward calculation shows that Eq.(4.3) becomes 



d 2 W a (z) 2{ 1 -\)dW a {z) X^- a > +1 



+ 



W a (z) = 0. 



dz 2 z dz z 

So, with v = ^2 > independent of 7 , we have 

d 2 W a (z) | 2(7-j) dWajz) W a (z) 
dz 2 z dz z a 

With another change of variable Y{z) = 7_1 W a (^) , (4.6) leads to 

d 2 Y 1 dY 



= 0. 



dz 2 z dz 



(7- I) 2 1 

z 2 z a 



Y = 0. 



(4.6) 



(4.7) 



Finally with the further change of variable £ = 2vz & , we have from (4.7) 

d 2 Y 1 dY 



+ 



de e da 



1 , PM7-1) 



r = o, 



(4.8) 



which is the equation of a modified Bessel function 35 of order 21/(7— 1) • The solution 
of Eq.(4.8) is 

W a (z)=a z 1_7 J 2 „( 7 -i)(2i/2 _ £) + c 2 z 1_7 i4: 2 ^ (7 _ 1) (2j/z _ ^), 

where / and if denote the modified Bessel functions of the first- and second- kind 
respectively 35 . The initial conditions lim W a (z) = 0, and lim W a (z) = 1 yields ci = 

and c 2 = r(C(7-i)) b y means of 



#„(*) « ir»| l 



(4.9) 



as z approach 0. Therefore, we have 

21/M7-1) 



W Q (*) = 



V^ Q (x;A) 



r(2i/( 7 -i)) 

2^(7-1) 



r(2i/( 7 -i)) 



2:1 7 ^2^( 7 -i)(2^z ^ or more explicitly 

r (7-i) x i-7^ (7 _ i)(2 ^VAx-^). 



(4.10) 
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Finally, we have for the (un-normalized) wave function (4.1) that 

^ A) = r(M7-l)) A ' (7 ~ lV ~ 7X M7-i)(2^-^)<M*). (4.11) 
It is quite clear by means of Eq.(4.9) that ]imip(x;X) = ipi{x) as expected. Conse- 

A — »0 

quently, the normalization constant N x of ip(x;X) must satisfy lim N x = 1. Some 
properties of the function K v (z) are in order 35 . The physical importance 35 of the 
function K„{z) lies in the fact that it tends exponentially to zero as z — » oo . The 
function K v {z) is defined, for unrestricted values of v , by the equation 

where 

w = (|rE t!r - + i +1) - 

fe=0 

The apparent discrepancy with (4.12) is resolved by the identity T(u)T(l -v) = sin ^ 7r) . 
For integer values or zero of v in (4.12), it should be understood that K n (z) = 
lim K v {z), where in this case 



fc=0 

1 7 °° C— ) fc 

+ (-ro(o)"Ew fc+1 )+^+fc + i)} ! 



(4.14) 



2 V 2 ; <^ LY ^ y ^ v /J fc!(n + £:)!' 

fe=0 



while 



*„(*) = "(c + log(-)KoW + £ + 2 + 3+ -- - + ^)- ^ 4 - 15 ) 

Here c is Euler's constant c = .57721 56649 The following identity will be also in 

used 

z< ^r = ~ vKv ^ - zK -+^ z )- ( 4 - 16 ) 

5. Lowest-order asymptotic perturbation corrections for a > 27 

In this section, we apply Theorem 1 and the trial function developed in section 4 
in order to obtain the eigenvalue perturbation expansions for the Hamiltonian (1.2). 
We consider first the case of a > 27 which leads to the divergence of the first-order 
correction of the regular Rayleigh-Schrodinger series. In a purely theoretical approach, 
Greenlee 2 has shown that the asymptotic perturbation expansion should take the 
form 

E (X) = E + E 1 X 2v ^ (5.1) 

valid for 27 < a and 2^(7 - 1) < 1. We note, for consistency, that we have re-produced 
the expression of Greenlee using our own notation. Eq.(5.1) is in complete agreement 
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with our prediction a > 2j or 2,(7 - 1) < 1 for v = obtained by means of Theorem 
1, as we shall show in this section. In order to apply Theorem 1, we need first the 
normalization constant N x of the trial wave function ip(x;\) , namely Eq.(4.11). This 
can be found by means of the condition || ij)(x;X) || 2 = 1 which leads to expression 



_ 2 4;/M7-l) ^(7-1) 

Nx [r(2K7-i)] 2 Jo 



r 2(l-7) 



^(7-1) (2,VAx 2I) 



tpi(x) 2 dx. 



(5.2) 



Lemma 1: For the ground-state, i.e. i = , we have 



r( 7 )r(i + 2,(7-1)) 



+ ( A M 7 -i)) 



(5.3) 



where v = and a > 2 7 . 



Proof: We note, by using (4.9) in (5.2), that N 2 m 1 . To find the order of the error term, 
however, we use the identity (4.12) which leads to K u {z) = £ ^(§)" 1 ' - £i ^(§ Y + ■■■■ 
Therefore 



^rd-Mj-i))] 2 ^-^^-^,,,, 



4, 2 

(5.4) 

For the ground-state, we have from (2.1) that ip (x) = ^p^af^ie^ 2 / 2 . Thus on 
substituting (5.4) into (5.2) we have, after some calculations, 

r ( 7)r ( i + 2,(7-1)) A +0{ } 

and the proof of the lemma follows by a very similar argument to that for Taylor's 
expansion of 1/(1-6). □ 
The reason of quoting the expansion (5.3) only up to order A 4l/(7_1) was guided 
by the error term in (3.12), as the following lemma indicates. 

Lemma 2: For the ground-state energy of the Hamiltonian (1.2), where a > 2 7 (or 
2,(7 - 1) < 1) , we have 



E (A) = 2 7 + 



16 z,M7-l);\2K7-l) + 1I 

r( 7 ) [r(2K 7 -i)P 



J x 1 -^ e- x " K 2v(l _ 1) {2vV\x-^)K 1 _ 2v(l _ 1) {2uV\x-^)dx + C^A 4 ^- 1 )). 



(5.5) 



Proof: Eq. (4.10) with (4.16) leads to 



dW a (x:X) 2v 2v{ -!-±) 



dx r(2,( 7 -l)) 



A K7-i)+^-^-7 Xi (2in /Ax-A). 



(5.6) 
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Furthermore, using Vo0*0 = yr^y^ 7 ^ e 2:2/2 , we find 



<W a (x;A) 
dx 



(7-|) 



cix 



2i/2^(7-i)A"(7-i)+i / 2 
x A r 1 _ 2l/ ( 7 _i)(2i/\/Aa; _ ^), 



a; 2 2^ e 



72 



(5.7) 



which leads to 



2(W a (x-\)Mx) 



dW a (x;X) 
dx 



(7-|) 



_d 



V'o(z)) = 



16 i/M7-l),\M7-l)+3 



r( 7 ) [r(2^( 7 - 

oo 

y x 1_ ^ e - K2 x 2y(7 _ 1) (2^VAx _ ^)i^ 1 _2^(7-i)(2^VAx _ 2T)rf a 



(5.E 



In order to use theorem 1, however, the trial wave function must be normalized. This 
is equivalent to multiplying (5.8) by the normalization constant N%, as given by (5.3). 
Now, since iV| is of order A 21 ^ 7-1 ) , out of the second term in (5.3) the multiplication 
allows us to have (5.8) as quoted, plus an error term of order A 4l ^ 7_1) as result of 
using (4.9). What remains is to show that the expression || [H + XV - E(\)]<j>\ | in 
(3.12) is also of order \ 2 ^~ l ) . This follows from (5.7) as follows 

dW a (x;X) f( 7-|) d~ 
x dx 



[H + \V- E ]tPo II = 2 



dx 



M*) II 



= 2 



8(/ 4,(7-l) A 2,(7-l)+l 



OO 

J x 1 ~-e~ x2 [K 1 _ 2v{l _ 1) (2uVX: 



x ^)] 2 dx 



1/2 



= 0{X 2v ^ 1 ~ r> 



where we have used Jr 1 _ 2l/(7 _ 1) (2iVAar&) w W-Mn-V) („^\ x -^)Mt-i)-i . The proof of 
the lemma then follows by use of Theorem 1, Eq. (3.12), and the variational estimate 
of E$ by means of (2.2). □ 

Because of the error term in (5.5), it is not necessary to compute the integral 
in (5.5) exactly but it is sufficient to estimate the integral using the asymp- 
totic series expansions of the modified Bessel functions #21/(7-1) (2^\/Ax~2t) and 
#i-2j/( 7 -i)(2^v / Ax - ^) by means of (4.12), up to the order cited. Since the order 
of the error term in (5.5) is A 4l ^ 7_1 ) while the integral is of order A 21 ^ 7-1 ) , we 
may consider, for fixed a, two regions < 2^(7 - 1) < \ and \ < 21/(7 - 1) < 1, or 
equivalently < 41/(7 - 1) < 1 and 1 < 4^(7 - 1) < 2 . For the first region, we have for 
the ground-state energy of the Hamiltonian (1.2) 

4(7 - l)i/ 4,/ (7-i)r(l - 21/(7 - 1)) 



E {\) = 2 7 - 



^2i/( 7 -l) + (3(^41/(7-1) j 



(5.9) 



r( 7 )r(i + 21/(7 - 1)) 

which follows from (5.5) using the asymptotic expansions of K 2v ^-i){2v\/Xx~'^) and 
#i-2 I >(7-i)(2^v / Ax - ^) by means of (4.9). In the case 7 = 3/2 (i.e. A = or a > 4), Eq. 
(5.9) yields 



E (X) = 3 + 



v^T(l + f) 



y + o{\ 2v ) 1 



(5.10) 
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where v = , as shown earlier by Harrell for the spiked harmonic oscillator Hamilto- 
nian (1.1). Important conclusion follows from (5.9). For a = 2(27- 1) or 21/(7 — 1) = 3 > 
we have 

E (\) = 2 1+ -^-yX + O(X). (5.11) 
r (7) 

This provides a single ground-state approximation formula for a wide class of Hamilto- 
nians H = -■£^+x 2 + ^ + ^, where a and 7 = l + ^y/i + AA are related by a = 2(27-1) . 
For example, for A = 0,i.e. 7 = 3/2, which yields a = 4, we have 

E (X) = 3 + -^VA + 0(A), 

as noted by Harrell. If a = 6, which implies 7 = 2 or A = 0.75 , we have 

£ (A) =4 + 2VA + 0(A). 

For the second region 1 < 41/(7 - 1) < 2 , or 27 < a < 2(27 - 1) , by using (4.12), we 
can easily show that 

^ (7 _ 1)) (2^--)^ 1 _ 2 , (7 _ 1)) (2^-- ) = r(2K7-i))r(i_- 2,(7-1)) ^ 

_ [T(2l/( 7 - l))] 2 l-4.( 7 -l) A l-2„(7-l) x ~h +2(7-1) + 

4(1-2,(7-1)) 
Consequently, (5.5) yields, for i < 21/(7 - 1) < 1 , 

E (A) - 2 7 + 4(7-1)^-^(1-2,(7-1)) ^-!) _ 2„r( 7 -£) 0(aM7 _ 1}) 

0( ) 7+ r( 7 )r(i + 2,(7-1)) (1-2,(7-1)) r( 7 ) A + U(A j ' 

(5.12) 

Again the result of Harrell for the Hamiltonian (1.1) follows for the case of 7 = | , i.e. 
A = , where, in this case, 3 < a < 4 or \ < v < 1 , and 

F m _o 1 4, 2 T(1-,) i^Ml A + OfA 2 n f513^ 

So(A) - 3+ ^m + u) A _ (i-,)v^ ( } " ( } 

For the rest of this section, we consider the case of 2,(7 - 1) = 1 . For this specific 
value Eqs. (4.10) and (5.6) read, for z = \~ v x , 

wm = ^ljr l( ^l) (5.14) 

^ - ^^(^), (5.15) 
dz 7 — 1 7 — 1 

respectively. Using the asymptotic expansions 

K (z) = [-c + log(2)-log(z)] + 0(z 2 ), #i(*) = -+0(*) (5.16) 

z 

which follow by means of (4.15) and (4.14) respectively, we obtain 



Perturbation expansions for a class of singular potentials page 14 



Lemma 3: For the ground state energy of the Hamiltonian (1.2), where a = 27 (or 
2^(7 - 1) = 1) , we have 



^o(A) = 27- (7 _ 1 1 )r(7) Alog(A) + 



-c(l + 7) + 21og(2( 7 -l)) 



A + 0(A 2 log 2 (A)), (5.17) 



( 7 - i)r( 7 ) 

where c = .57721 56649 ... is Euler's constant. 
Proof: We should note first, in this case, 

|| [H + XV- £(A)]Vo H 2 = 0(A 2 log 2 ( 7 )) 

which follows from 

o\2 poo /T 

|| [H + XV - 25(A)]*, || 2 = (7 _ 1)2r(7) I x^e-* 2 Kt(^- iX ^)d X 

by use of (5.7). Since we are only interested in finding the order in terms of the 
parameter A , the problem reduces to a search among the smallest value of A 2 log 2 (A) , 
-A 2 log(A), and A 2 , for small values of the parameter A. Therefore for sufficiently 
small A we have || [H + XV - E(\)]ip || 2 = 0(A 2 log 2 (A)) as noted . What remains is to 
compute 



„/ TI r , , x , / n dW a (x; A) 
2(W a (x;X)Mx), 



(7-j) d 
x dx 



4A3/2 



( 7 -i) 2 r( 7 ) 



Jo 7-1 7-1 



by using the asymptotic expansions (5.16) up to the order A 2 log 2 (A). The lemma 
then follows after some straightforward calculations. It is important to note that the 
normalization constant A^a, as given by (5.2), yields in this case 

N- 2 = 1 + |2^ 2 (-1 + 2c + 21ogH) + ^cJjA- + ^AlogA + . . . 

and will contribute to the error term in a similar manner to that mentioned in lemma 
2. □ 

The results of Harrell, the case a = 3, follows immediately from (5.17) for the 
case of A = (or 7 = 3/2 ), i.e. v = 1 , namely 

E (X) =3-^Alog(A)-^A + 0(A 2 log 2 (A)). (5.18) 

It is clear that these expressions are valid for A much smaller than unity. 

6. Lowest-order asymptotic perturbation corrections for 27 > a > 7 + 1 

In this section, we discuss the case of 27 > a > 7 + 1 or equivalently the case 
of 1 < 2^(7 — 1) < 2 . It is clear by now that, for 1 < 2^(7 - 1) < 2 , the first-order 
Rayleigh-Schrodinger corrections exist but the second-order corrections diverge. Thus, 
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the improved perturbation procedure gives an explicit term between the first and the 
second order. Let us first consider the case of 1 < 21/(7 _ 1) < 2 , here we rely on the 
asymptotic expansion of the modified Bessel function K as given by (4.12). We note 
from (4.5) and (5.6) that 



[H + XV- £(A)]Vo 



'■)dx 



32^(7-1) A 2-(7-l)+l ,00 

r( 7 )[r(2,( 7 -i))P 7 x ^ ^(7-d-i(^ - 
as a consequence of the known identity ^(2) = K_ v {z) . Using (4.9) we have 

|| [ffo + AF-S(A)]^ || 2 -0(A 2 ) (K2i/(7-l)<2), (6.1) 
which leads to the following lemma. 

Lemma 4: For the ground-state energy of the Hamiltonian (1.2) where a > 7 + 1 (i.e. 
1 < 21/(7 - 1) < 2) , we have 

16 j/M7-l),\M7-l) + 2 



E (A) = 2 7 + 



r( 7 ) [T(2i/(7-l)] s 



y x 1 "^e" x2 ^ (7 _i ) (2^v^x"^)ii'2^ (7 _i ) _i(2^v^\a; _ *)dx + 0(A 2 ). 



(6.2) 



The proof of this lemma is similar to that of lemma 3, therefore we omit it. The 
computation of the integral in (6.2) up to the order of A 2 yields the perturbation 
expansion 



0{ ) 7+ ^(7)^(7-1)) 



(21/(7 -i)-i) r( 7 ) 



A + 0(A 2 ), (6.3) 



as the result of 



1 , r(2i/( 7 -i))r(i- 21/(7-1)) 



j^VAx-*)- 2 ^- 1 



(isy/Xx 2I ) 2 
T(2- 21/(7-1) 



-(iVA^) 2 "^ 



r(i - 21/(7 - 1)) 

(iz-s/Aa; - ^) 2 



r(i + 21/(7 - 1)) r(2 + 21/(7-1) 



The result of Harrell 1 follows immediately from (6.3) in the special case A = 0, 
namely 

(6.4) 



£ (A) = 3 + 



4 " 2 " r ' 1 -"»A- + 2nM)A + 0,A-), 



V5rT(l + f) 
where 1/ = ^ an d 5/2 < a < 3 . 

For the case 2^(7 - 1) = 2 , the norm || [H + XV - E(X)]ip || can be computed easily 
by means of Eq.(4.5), which yields 



[Ho + XV- E ]^o 



dW a {x;X) 



dx 

\3 



(7- \) d 

x dx 



2VA 



(7 



^^/^e^[^(^-)]^ = 0(A 2 ) 
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as a consequence of 

Ki(z) = \ + + 2c) + i(- log(2) + log(z)))z + 0(z 2 ). 

Furthermore, since the trial wave function takes the form 

ip(x; A) = 2v 2 \x~i K 2 (2uVXx~i : )tpo{x), 
we may compute the normalization constant Eq.(5.2) as 

^ = l + ^A+^A 2 log(A) + 0(A 2 ). 

in a similar fashion to the proof of lemma 1. Harrell, in his investigation, claims that 
the ground-state eigenvalue perturbation expansion for the spiked harmonic oscillator 
Hamiltonian (1.1), a = 5/2, is given by 

E (X) = 3 + ^$^X + ^=X 2 log(A) + 0(A 2 ). (6.5) 
V 71 " v 71 

The following result does not confirm his claim, but shows that it is slightly different 
even by means of Harrell's own methodology. As our calculation will show, the ground- 
state perturbation expansion in the case of a = 5/2 is actually given by 

E (X) = 3 + —^-X + —j=X 2 log(A) + 0(A 2 ), (6.6) 
V 71 " V 7T 

with a multiple of 2 in the log term in contrast with (6.5). In order to verify (6.6), 
we adopt two different approaches, first we use Harrell's method then we apply our 
generalization. We have for a = 5/2 , the asymptotic perturbation expansion according 
to Harrell reads 

E (X) = 3 + 2(W 5/2 %, - ^M>) + 0(A 2 ), 

where ip (x) = -^jjxe^ x2/2 . For W 5/2 , Harrell used the asymptotic approximation of 
the modified Bessel function 

to show that 

W 5/2 (z) = 1 - 4= + - 4(4(c + log(2)) - S)z + . . . for large z. 

v 2 z 

If we differentiate W 5/2 (z) with respect to x, keeping in mind z = X~ 2 x, we have, 
after some calculations, 

oo oo 

om, / dw 5/^A d. f . 16A f 3/2 2 64A 2 logA f _2 
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which yields (6.6), since J xe ^ dx = \ and / x 3/2 e 2,2 dx = |r(|) . Numerically, how- 

o o 

ever, Eq.(6.5) is more appealing than (6.6) for a wider range of the parameter A 
smaller than unity since (6.6) reduces the applicable range of A by almost one-half. 

Lemma 5: For the ground-state energy of the Hamiltonian (1.2), where 21/(7- 1) = 2 , 
we have 

E (X) = 2 7 + ^A + ^A 2 log(A) + 0(A 2 ). (6.7) 



Proof: For 21/(7- 1) = 2, we have, using Eq.(6.2), that 

_ 00 

E (X) =2 7 + ^ / a; 1 -ste- a;2 X 2 (2z/v / Ax-*)^i(2 l /VAx-*)da ; + 0(A 2 ). 

1 (7) J 



By means of the asymptotic expansions 

1 ,1 , _ > 1 . ,1 



and 



= § + (^(-1 + 2c) + 2 log(-z))z + 0(z 2 ) 



^) = ^-5 + (^(§-2c)-ilog(^ + 0(A 



where z = 2v\f\x & , we have, after some calculations up to order A 2 , that 

00 00 
E (X) = 2 7 + y a^e-^dz + i^A 2 log(A) | xef'dx + 0(X 2 ). 


oc 

This leads to (6.7), since / x 1+ ^er x2 dx = ±r(l + . □ 



7. Further Cases 

The expansions developed above can be extended to the case of 21/(7 - 1) > 2 and 
the region of a < § which was not studied by Harrell 1 . For example, in the case of 
2 < 21/(7 — 1) < 3 , the second order of the perturbation correction exists but the third 
order diverges. By using the modified perturbation theory we find 

W-27 + r(7)r(1 + 2l/(7 _ 1)) + (2K7-D-D r (7 ) A 

2(3 - 4i/( 7 - i))j/ 3 r( 7 - i) , , , N 

v - - - v - A + 0(A 3 ). 



(2-2K7-l))(l-2K 7 -l)) 2 r(7)' 
Again in the case of 7 = 3/2, HarreH's formula should read, for 7/3 < a. < 5/2 , 

w) = 3 + + r r(i :^A - 4 ; 3 -y r( t^A 2 + tw, (7. 2) 

V7 V^T(l + ^) (i/-l)V5F (2- -z/) 2 V^ K h K ' 
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where v = . Furthermore, for 3 < 2,(7 - 1) < 4 , we have 

F (\\ -*v + 4(7-1)^(7-1)^1-2,(7-1)) ^-!) , 2^(7-^) 

^ o(A) - 27 + r( 7 )r(i + 2K7-D) A + (2^(7 -D-D r( 7 ) A 

2(3- 4,(7 -i)), 3 r(7-i) _ a2 2^(7-^) A 3 



(2 - 2,(7 - i))(i - 2,(7 - i)) 2 r( 7 ) r( 7 )(i - 2,(7 - D) 2 (2 - 2,(7 - 1)) 

+ 0(A 4 ). 

(7.3) 

Similar expressions can be obtained for n < 2,(7 - 1) < n + 1, n = 4, 5, . . . , as well for 
2,(7 - 1) = n. It is important, however, to note that there are an infinite number of 
cases as n increase to infinity for a > 2 . For example, if we restrict the value of 7 
to 3/2 , then in this case 2 < a < 5/2, and obtain an infinite number of perturbation 
expansions consisting of analytic parts of degree A™ in addition to one correction term 
in each expansion as 2 + < a < 2 + \ . Similarly, the cases a = 2 + n = 3, 4, . . . , 
need special treatment, as we have mentioned above. It is interesting to note that the 
A-term in (7.1) or in (7.3) is identical with the corresponding A-term in the Rayleigh- 
Schrodinger series. This can be easily verified by a comparison of the coefficient of A 
in (7.1) or in (7.3) with the coefficient of the A-term in (2.3). However this is clearly 
not the case for the A 2 -term. 

8. Excited states 

The perturbation expansions developed so far were restricted to the ground-state 
energy, however, it is a matter of calculation to extend these results to the excited- 
state energies. First, it should be noted that the exact solution of the unperturbed 
part of the Hamiltonian (1.2) has very little to do with the order of the error terms 
in the perturbation expansions. Therefore it is expected that the order of the error 

terms remains the same for the excited states, i.e. i = 1, 2, Following the discussion 

in Sec. 3, the asymptotic expansions for the eigenvalues for excited states are given 
by means of Theorem 1 as 

"(7-|) d 



Ei{\) = 2(2* + 7) + 2(W a (x; \)^{x), d ' 1 " ( ''' : ' ' 



dx 



MX)) (8.1) 



dx 

+ 0(|| [H + XV- E*]W a (x; A)^ || 2 ), 

where i[>i(x), i = 0,1,2,... are given by (2.1) and the energy of the unperturbed 
Hamiltonian Ef is given by (2.2). In order to compute (8.1) explicitly, we notice first 
that 



dW a {x;X) 
dx 



(7-|) d 
x dx 



... 16( 7 ) t ,M7-i) A 2K7-i)+| 

^ [x)) ~ iir( 7 ) [1W7-DP 



x / x 1 2 »e 



o 



A + l) 7; s 2 )iFiH; 7 + D^ 2 ) - [iF^-i; 7; x 2 



7 

x K 2u ( 1 -i)(2v^x~^)K 1 _ 2v ( 1 -\)(2v\/\~x~^),dx 



(8.2) 
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where iF 1 is the confluent hypergeometric function mentioned earlier Eq.(2.4). In 
producing (8.2) we have used the following identity 35 

^iFi (a; b; z) = ^(a; b; z) - ^Z^. lFl (a; b + 1; z). 
For the case of < 2^( 7 - 1) < § or a > 2(2 7 - 1) , we find, by means of (4.9), that 



2(W a {x-\)^{x) 



dW a (x;X) 



dx 



x j xe 




x dx 
2(7 + *) 



. , .,_ 4( 7 ), r(l- 2i/(7-l)) 

4l/(7-1) \ 2i/(7-1) 

^ iW; ~i!r( 7 ) i/r(2i/( 7 - 1)) 



7 



*!(-*; 7; x^i^-i; 7+I; x 2 ) - (-i; 7; x 2 )]< 



(8.3) 



This result required the investigation of some integrals of the type 

00 00 

J xe~ x2 iFi(-i]j;x 2 )iFi(-i]~f + l;x 2 )dx &nd J xe - * 2 [ii<i(-z; 7; x 2 )] 2 dx. 



Lemma 6: For d > and s > 

00 

/ ^e-^iF^asfc; kt^F^a' ]b' ;k't)dt = S - d T{d) 



m=0 



Proof: From the series representation, Eq.(2.4), of the confluent hypergeometric series 
1F1 , namely 

1 F 1 (o;6;fct)=f;^^ and lFl (a>; b'; k't) = £ 

^„ (o)m m! ^ (&')„ n! 

m=0 v ' n=0 v ' 



we have 



00 

/oc 00 
t d - 1 e- s \F 1 (a;b;kt) 1 F 1 (a';b';k't)dt= ^ ^ 



m— n— 

00 00 



e - s t t d+n+m-l dt 



EE 



(a) m (o') n k m k' A )n{ l r{d + m)nT{d + m) 



-- d E E 



m=0 



(a')n(rf + m) n {—) 



n=0 



(b')r. 



ni 



(«Uf) r 

(b)m 



■T(d + m) 



-d-ni j\ 771 / / j 1 ul ^' \ ( a )m(d) m ( s ) 



m=0 



with > and s > , and where we have used the series representation of 2 F\ as 
given by Eq(2.4) and the Pochhammer's identity r(d + m + n) = (d + m) n T(d + m) = 
(d + m) n (d) m T(d) . □ 
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As a consequence of this lemma, we have 

oo 

J ze^ViH; 7; s 2 )iFi(-i; 7 + 1; x 2 )dx = ^i^- 3 f 2 (-t, 1 - 7, 1 ;7 , 1 - 7 - i; 1) 



and 



00 

J xe~ x2 [rFiH; 7; x 2 )] 2 dx = ilz2ll 3 F 2 H, -7, 1; 7 , 2 - 7 - i; 1), 



(8.4) 



(8.5) 



where the Chu-Vandermonda identity 



2 Fi(-n,a;c; 1) 



(c - q) T 

(C)n 



has been used. Therefore, from (8.3), we have 



dW a (x; A) 



(7-i) 



da; 



Ml = 4 (7)» r(l- 2^(7-1)) 4, (7 -i) A 2,( 7 -i) 



t F 2 (-i, 1 - 7,1;7,1 - 7- i; 1) 



(7-1), 
2(7), 



iir( 7 ) i/r(2i/( 7 - 1)) 
3^2(-i,-7,i;7, 2 -7-i; 1 ) 



i.6) 



For < 2i/( 7 - 1) < \ or a. > 2(2 7 - 1) , we have 

E (X) - 2(2i + >y) + 4(7)i r ( 1 - 2 ^- 1 )) ,^-D;2>h-i) 
^(Aj-A^ + 7)+ ., r(7) ^ r(2iy(7 _ 1)} " 

(7-l)i 

3 F 2 (-i,l - 7,1; 7, 1-7-2;!) — —^ 3 F 2 (-i,-7,l; 7,2-7 -i;l) 



2(7), 



(8.7) 



+ 0(\ 



4 1/(7-1) 



), i = 0,l,2,.. 



Similar expressions can be obtained for the other cases by means of lemma 6. An im- 
mediate extension of Harrell's expansions to excited states i = 1, 2, . . . for the Hamil- 
tonian (1.1) can be obtained by setting 7 = | in (8.7). 

9. Conclusions 

In this article, we have applied and extended Harrell's modified perturbation the- 
ory to treat a wider class of singular Hamiltonians given by (1.2). Our extensions 
allow us to recover Harrell's formulas for the spiked harmonic oscillator Hamilto- 
nian (1.1) as special cases. Further, we were able to extend Harrell's results to the 
exited-state energies, again as special cases of our general treatment. We have also 
now corrected the perturbation expansion (6.5) for the case a = § : this formula has 
been used without correction since the very early work of Harrell. Some interesting 
questions which remain to be answered are as follows. Is the modified wave function 
(4.11) sufficient to extend the perturbation expansions presented here to higher or- 
ders, or will further modifications need to be introduced? Why can the second-order 
corrections in (7.1), (7.2), and (7.3) not be recovered from the corresponding terms 
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in the regular Rayleigh-Schrodinger series? Does this fact indicate that the trial wave 
function indeed requires additional modification? We hope that the present work will 
encourage further research into this interesting class of singular Hamiltonians. 

Acknowledgment 

Partial financial support of this work under Grant Nos. GP3438 and GP249507 
from the Natural Sciences and Engineering Research Council of Canada is gratefully 
acknowledged by two of us (respectively [RLH] and [NS]). 



Perturbation expansions for a class of singular potentials 



page 22 



References 

1 E. M. Harrell, Ann. Phys. 105, 379 (1977). 

2 W. M. Greenlee, Bull. Amer. Math. Soc. 82, 341 (1976). 

3 R. Hall, N. Saad and A. von Keviczky, J. Math. Phys. 39, 6345 (1998). 

4 R. Hall and N. Saad, J. Phys. A: Math. Gen. 33, 569 (2000). 

5 R. Hall and N. Saad, J. Phys. A: Math. Gen. 33, 5531 (2000). 

6 R. Hall, N. Saad and A. von Keviczky, J. Phys. A: Math. Gen. 34, 1169 (2001). 

7 R. Hall, N. Saad and A. von Keviczky, J. Math. Phys. 43, 94 (2002). 

8 Omar Mustafa and Maen Odeh, J. Phys. B: At. Mol. Opt. Phys. 32, 3055 (1999) 

9 Omar Mustafa and Maen Odeh, J. Phys. A: Math. Gen. 33, 5207 (2000). 

10 J. Skibihski, e-print |quant-ph/0007059| . 

11 J. R. Klauder, Acta Phys. Austriaca Suppl. XI, 341 (1973). 

12 J. R. Klauder, Phys. lett. B 47, 523 (1973). 

13 J. R. Klauder, Science 199, 735 (1978). 

14 H. Ezawa, J. R. Klauder, and L. A. Shepp, J. Math. Phys. 16, 783 (1975). 

15 B. Simon, J. Functional Analysis 14, 295 (1973). 

16 B. DeFacio and C. L. Hammer, J. Math. Phys. 15, 1071 (1974). 

17 T. Kato, Prog. Theor. Phys. 4, 514 (1949). 

18 T. Kato, Prog. Theor. Phys. 5, 95 (1950). 

19 T. Kato, Prog. Theor. Phys. 5, 207 (1950). 

20 L. C. Detwiler and J. R. Klauder, Phys. Rev. D 11, 1436 (1975). 

21 Yudell L. Luke, The Special Functions and their Approximation, Vol. I (Academic 
Press, 1969). pp. 43 

22 T. Kato, J. Phys. Soc. Japan 4, 334 (1949). 

23 Evans M. Harrell II, Proc. Amerc. Math. Society 69, 271 (1978). Theorem 2 

24 G. Temple, Proc. Lond. Math. Soc. 29, 257 (1928). 

25 G. Temple, Proc. Roy. Soc. 119, 276 (1928). 

26 D. T. Hayes, Canadian J. Phys. 49, 218 (1971). 



Perturbation expansions for a class of singular potentials 



page 23 



27 M. Znojil, Physics Letters A 164, 138 (1992). 

28 M. Znojil, J. Math. Phys. 30, 23 (1989). 

29 V. C. Aguilera-Navarro and E. Ley Koo, Int. J. Theor. Phys. 36, 157 (1997). 

30 P. Chang and Chen-Shiung Hsue, Phys. Rev. A 49, 4448 (1994). 

31 V. C. Aguilera-Navarro, A. L. Coelho and N. Uttah, Phys. Rev. A 49, 1477 (1994). 

32 R. Hall and N. Saad, Canad. J. Phys. 73, 493 (1995). 

33 V. C. Aguilera-Navarro, Francisco M. Fernandez, R. Guardiola and J. Ros, J. 
Phys. A: Math. Gen. 25, 6379 (1992). 

34 M. J. Jamieson, J. Phys. B: Mol. Phys. 16, L391 (1983). 

35 M. Abramowitz and I. A. Stegan, Handbook of mathematical functions (Dover 
Publications, Inc., New York, 9 th Edition (1970)). pp. 507, formula 13.4.12 



